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Abstract 
00 _ 

^D I Neutrino fluxes at telescopes depend on both initial fluxes out of astronomical bursts and flavor 

^N ' mixing during their travel to the earth. However, since the information on the initial composition 

1^ ■ requires better precision in mixing angles and vice versa, the neutrino detection at telescopes 



for itself cannot provide solutions to the both problems. Thus, a probability to be measured at 
long baseline oscillation is considered as a complement to the telescope, and problems like source 



r^ ■ identification and parameter degeneracy are examined under a few assumptions. 
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I. INTRODUCTION 

Cosmic neutrinos can be classified into stellar, galactic, and extra-galactic neutrinos ac- 
cording to astronomical source. Stellar neutrinos include solar and supernova neutrinos, of 
which energy scale is order of (l-lO)MeV since they are mainly produced by nuclear inter- 
action. There are fair records in which they were identified as from the sun and SN1987A 
and examined for oscillation and matter effect [l||2]- Ultrahigh energy (UHE; > lO^'^eV) 
cosmic particles are regarded to have their origin in extra-galactic source like active galactic 
nuclei(AGN) and some of gamma ray bursts(GRB). Mechanisms to accelerate protons to 
high energy have been searched, from (9(PeV) to GZK limit, in GRB |3|]|j]. In astronomy, 
GRB's are detected once a day on average. If an accelerated proton produces a pion, about 
20% of its energy is transferred to neutrinos, although more than one pion can be produced 
from a proton. Thus there are abundant source of neutrinos whose energy is higher than 



Q. 



lOOTeV 

Decays of pions, such as vr^ -^ Z^"*" + t'/i -^ e"*" + z/^ + i/g + ^^i or vr" ^ yU~ + z/^ — ^ 
e~ + z/^ + z/g + t'^, are the main process of neutrino production. The initial neutrino flavor 
ratio can be determined, depending on the charges of pions and the energy loss rates of 
muons in processes of pion decays. Whether pions are produced by pp or by py collision, 
the charges of pions and the ratio of z/g/z/g are different. Although both pions and muons 
lose energy in environment that pions were borne, long-lived muons are more probable to 
interact with environment before decaying. The pions that decay into neutrinos without 
muon decays are called "muon-damped source" . Such electromagnetic energy loss of muons 
for GRB becomes significant when the energy of 7-ray E^ is > lOOTeV and, thus, the 
energy of neutrino E^, produced from muon-damped source is > ITeV 6|. The intensities of 
neutrinos in lO^^eV < E < lO^^eV and higher energy were studied. Its model-independent 



a 



upper bound was predicted |^[3l. If the decay mode of a pion includes muon decay, that 
will be called simply "pion source" in comparison with muon-damped source. The initial 
neutrino composition $°(z^e) : ^^{'^ij.) '■ ^^{'^t) = 1:2:0 from the pion source is gradually 
replaced by the composition $°(z/e) : $°(t'^) : ^^{lyr) = 0:1:0 from the muon-damped 
source when the energy increases passing about ITeV 6||8|. When the main source of 
neutrino production in atmosphere is pion decay, the shift of the ratio $°(t'^)/$'^(z/e) from 
2 to infinity according to increasing energy [91] [10| was tested in SK [ll| [12] . The transition 



in the flavor ratio can be naturally assumed to occur at the region of GRB's, and the 
corresponding aspect was discussed in Ref . [6| . The relative flavor ratios of detected fluxes 
vs. energy exhibit saturation and transition. The saturation of <l>*(t'^)/<l>*(z/e) close to one 
implies that the beam originated from pion source for lower energy, while the transition of 
$*(z/^)/$*(z/e) to a larger one implies that the beam includes portion of muon-damped source 
for higher energy. 

There are a few neutrino telescopes under construction, which can detect high energy 
neutrinos oi E,y > O.lTeV, e.g., IceCube [l3|, Antares pj|, etc. They are designed to detect 
VHE and UHE neutrinos and to distinguish a flavor from others. With a future Cerenkov 
detector, u^ is easy to be detected due to the long tracks of muons, while u^- is distinguishable 
from z/g by double-bang event only near PeV range. At (9(PeV), the neutrino telescopes will 
be able to observe both double showers of z/,-, one due to r production and the other due to 
r decay, and W~ resonant event at 6.3PeV that identifies Ue |15l |. 

For cosmic neutrinos, the oscillation factor in a transition probability is averaged out 
due to long distance of travel and high frequency. Thus, the telescope experiment can be 
effective measurements of mixing parameters, provided that the neutrino flavor ratio at the 
source is known [l6| The current data of mixing parameters are phrased by the ranges in 
the magnitude of PMNS elements; 



If/, 



PMNS] 



'^ 0.79 -0.86 0.50-0.61 0-0.20 ^ 

0.25-0.53 0.47-0.73 0.56-0.79 
\ 0.21 - 0.51 0.42 - 0.69 0.61 - 0.83 / 



(1) 



at 3cr level jl7| . The matrix implies the following values of individual parameters; Am2i 
(7.1 - 8.9) X 10-5eV^ jAm^J = (2.2 - 3.0) x lO'^eV^ sin^ ^12 = 0.24 - 0.40, sin^ ^23 



0.34 — 0.68, and sin^ 6'i3 < 0.040, all at 3a level [l8|]. The type of the yet-undetermined 
among neutrino parameters can be organized case by case: First, normal hierarchy (NH) or 
inverse hierarchy (IH): the sign of Am'^i = m| — ml is unknown while the global best-flt is 
obtained as |Am|^| = 2.6 x 10^'^eV^. Second, the sign of ^23~7r/4, i.e., whether the majority 
of Ut- is the heavier(z/3 for NH) or the lighter(z/2 for NH) is unknown. Third, a certain value of 
oscillation probability has inflnite number of candidate combination of 6*13 and 5cp [l9| [20|] . 
The ambiguity due to the above degeneracies needs to be clearly distinguished from that 
due to uncertainties. 



In Sec. ini the probability of a neutrino oscillation is reviewed, focused on how the degen- 
eracies are generated and which oscillation experiment each mixing angle is most sensitive 
to. In Sec. Illlt examined are the sensitivities of neutrino fluxes to neutrino mixing angles. 
In Sec. IIV[ we discussed a few points including parameter degeneracies and source identifi- 
cation, where the data of neutrino telescopes are considered in company with the data of a 
long baseline (LBL) oscillation under a few assumptions. Concluding remarks follow in Sec. 

El 

II. REVIEW ON DEGENERATE PROBABILITY OF OSCILLATION 

The probability of neutrino oscillation may be said to be degenerate because different 
sets of parameters result in the same value. For example, the probability of transition from 
Ua to Ujs in two-neutrino oscillation with a single mixing angle and a single mass-squared 
difference Am'^ 

Pap = 5^p - {25^p - 1) sin^ 20sin2(^^), (2) 

is invariant under switching the sign of Am^ or changing the angle 9 with its complementarity 
angle tt/2 — 9. Since the two sets, (6', Am?) and {tx/2 — 9^ — Am^), are not physically different, 
the above probability is twofold degenerate simply due to {9, Am?) and (7r/2 — 9, Am^) or 
due to {9, Am^) and (6*, -Am^). 

The oscillation probability extended to three neutrinos in vacuum 

J2 R4Uc.^u;,u:^Up,] sin' i^^) 

J^ A AmlL 

± 2^ 5] Im[U^,U*^,U:^Uf3,]smi^^) (3) 

i—l j—i+1 

is given in terms of 3 x 3 unitary transformation matrix U and three mass-squared differences 
Ar7i|j = m'j — r/i?, where each of v^ and v^ may be one of z/g, i/^ or v^. In effect Ar7i|^ ~ 
Am32 so that they will not be distinguished hereafter. The PMNS matrix in standard 
parametrization is given by 

C12C13 S12C13 siae 



C/a, 



-S12C23 - Ci2S23Sl3e*^ C12C23 - 3125235136"^ S23C13 
S12S23 - Cl2C23Sl3e*^ -C12S23 - Sl2C23Sl3e"^ C23C13 



where Sij and Qj denote sin 6ij and cos 6ij with the mixing angle 6ij between i-th and j'-th 
generations, respectively, and 6 denotes a Dirac phase. 

The degeneracy of probability could be stemmed from three mixing angles and three 
mass-squared differences. However, since a number of successful experiments found the 
values of some physical parameters, the possible multiplicity of degeneracy at present is at 
most eight. A twofold degeneracy in the eightfold degeneracy has its origin in the ambiguity 
of a value of 623 from its complementary angle, tt/2 — 623. The z/^-disappearance probability 



which is most sensitive to determine the value of 6*23 is expressed by 20 1 

1 - P{u^ -^ lyf,) = sin2 2023 cos2 ^13 sm'( '^^^^ ) 



4E,, 



sin(— J^) sin(— J^) sin^ 2^23 



4^^ ' ' AE^ 
■ (cos^ 013 cos^ 012 — sin 0i3 sin^ 023 sin 20i2 cos S) , (4) 

where the leading two terms are invariant under the exchange of 623 with its complementary 
angle. On the other hand, the i/g- appearance probability with matter effect A, 

P{v^{V,) -> v,{V,)) = sm^ 023 sm^ 20i3 \\ ^ ^'/^^y 2 

2 



+ - — ^ sin 2023 sin 20i3 sin 20i2 cos((5 ± ATO31 LlAE^) 
Amg]^ 

sin((ATO§i =F A)L/AE) sai{AL/AEy) 



(1TA/Amij(^/Am2,) 

.^-5-) cos 023 Sm 2012 (^/^^2^). 



+ (3lfcos^023sin^20i2 7-;;:r" , (5) 



provides multi solutions to (6^13, 5) pair and ( Am^ i , 5) pa ir, where A = 2\/2G-pYcpEy is 



given with a density p and an electron fraction Y^. 19|] 20|]. The leading term in Eq. ([5]) 
is not sensitive to the sign of Ar/ig^, if the matter effect A is not significant. The second 
term is also not indicative of the sign due to cos 5 and a suppressing factor ^in\^/ ^m\^. So 
the (Ar?7,3]^, 5) pair causes another double degeneracy to the probability. It is well known as 
a matter of normal hierarchy or inverse hierarchy. The multiple possibilities in the pair of 
(6'i3, 5) for a value of sin 29i3 cos 5 make up the eight-fold degeneracy with other two double 
degeneracies. 

If astronomical neutrinos from origins like grb's or other types of extragalactic bursts are 
considered for the flavor transition, in the limit L ^ 00, the probability in Eq. ([3]) reduces 



to 

2 3 

P*^ ^ 5^p~2j2Yl Re[U^.U;,U:^Up,] (6) 

= x^if/..rit//3.r. (7) 

Since the probability of relatively long-distance oscillation compared to wavelength is av- 
eraged out, the dependency in the sign of Am^j is hidden, so that the probability is blind 
to Am^j. Since it depends on only the absolute values of the elements in PMNS as shown 
in Eq. ([7]), its dependency on the phase 6 is simply cos 5. Sensitive dependencies between 
other mixing angles and P^a will be discussed in detail in Section IV. 

Degeneracy problems are caused since the number of effective measurements is not suffi- 
cient to specify all the physical parameters and the forms of probabilities are multi- variable 
sinusoidal functions. Here, the ambiguities from three types of degeneracies, the sign of 
Am|^ 211], the sign of tt/A — 623 [22], and different pairs of (^i3,(5)'s [23] will be clearly 
distinguished from the ambiguities from uncertainties. The degeneracy is a problem caus- 
ing ambiguity even when an average probability or an average flux without uncertainty is 



applied to determine parameters. For more review, see references 24] 



III. LIMIT OF NEUTRINO TELESCOPE 

The initial flux of astronomical high energy neutrinos is assumed to be attributed from 
the decay of pion produced in py collision. The relative ratio of neutrino flavors from 
pion decay is ^^{lye) '■ '^'^(^^) : '^'^(^r) = 1:2:0. If the daughter muon in pion decay 
does not decay due to the electromagnetic energy loss, the flavor ratio is changed into 
$°(z/e) : ^^(i^fj) : $°(z^r) =0:1:0. So, the composition ratio $°(z/^)/$°(z/e) = 2 and its 
transition to inflnity as the energy of neutrinos increases may occur in region of GRB's, as 
does in atmosphere. However, it is impossible to estimate the initial flux of neutrinos as 



being produced in GRB's. One can but conjecture it from detected fluxes at telescopes [25] . 
If the neutrino mixing is of tri-bi maximal type as the simplest example, the pion source is 
identified by the fiux ratio at telescopes, $*(z/e) : '^*(t'/^) : "^^l^r) = 1:1:1, while the muon- 
damped source is identified by the fiux ratio at telescopes, 1 : 1.8 : 1.8. Thus, the transition 
in fiavor ratios of detected fiuxes at GRB's occurs with the relative ratio $*(z/^)/$*(z/e) from 1 



to 1.8 as the energy of neutrinos increases. In reality, it is impossible to determine, simply by 
reading results at telescopes, whether initial beams are pion source or muon-damped source, 
without any assumption, for the following reasons: First, any information from astronomical 
neutrino bursts on the composition in a mixture of two types of sources cannot be obtained 
without using telescopes. Second, the mixing angles and masses are still bearing too broad 
uncertainties to analyze the fluxes, and moreover the fluxes to be detected at telescopes are 
significantly sensitive to mixing angles. Hereafter, we took a strong assumption which is 
that pion source and muon-damped source can be distinguished, in order to illustrate the 
point of the sensitivity to mixing angles. 

Unless neutrinos decay, <l>°(z/e) + '^°(t'/x) + *&°('^r) = ^\i'e) + ^\i^fj.) + ^\i^T) and J2a^a(3 = 
1, where $*(z^/3) = '^a^a/s^'^i^a)- Since relative fluxes at a telescope can be normalized 
such that X]a'^°('^a) — Xla'^^l'^a) — ^5 ^^c normalized flux $*(z/q,) is a linear combination 
of -Pga's, for example, $*(t'Q,) = P^„ for muon-damped source and $*(t'Q,) = 1/3{P^^ + 2P^^) 
for pion source. The symmetric matrix P^^ in Eq. ([7]) has only three independent elements. 
Thus, when the initial condition of neutrino beam is assumed such that both pure pion source 
and pure muon-damped source are allowed, $*(//,-) and one of the following four fluxes may 
be unnecessary: $*(z^e) from pion source, ^^{i'^) from pion source, $*(z/e) from muon-damped 
source, and $*(t'/x) from muon-damped source. In reality, either $*(t'e) or $*(t'/x) is more 
likely to be a mixed flux from pion source and muon-damped source, rather than a flux from 
the pure pion source or the pure muon-damped source. In order to examine the sensitivities 
to mixing angles, however, we consider the pure pion or the pure muon-damped source. 
FIG. [Hand FIG. |2] illustrate the dependence of $*(t'e) on 6*23 and the dependence of $*(t'^) 
on 623, respectively. For a given ^13, the range of S from to 2tt makes the curves thick 
bands. The red (upper in Fig. 1 and lower in Fig. 2) and blue (lower in Fig. 2 and upper in 
Fig.l) bands indicate the fluxes from the pure pion source or from the pure muon-damped 
source, respectively. 

The symmetric probability in Eq. ([71) has ranges estimated from PMNS at 3a CL in Eq. 

(HD; 



^ 0.48 - 0.64 0.12 - 0.34 0.11 - 0.35 ^ 



Pa/B 



y/ 0.33-0.53 0.30-0.41 

\ V V 0.33- 0.47 y 

Thus, the fluxes accompany nontrivial ranges. The plots in FIG. [T]and FIG. [2] illustrate the 



(8) 




0.65 070 0.75 0.80 0.85 90 0.35 



923 



0.65 0.70 0.75 0.80 0.85 0.50 95 



FIG. 1: ^\i^e) vs. 6*23 for fixed values of 6*13. 
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FIG. 2: $*(fj^) vs. 6*23 for fixed values of 613. 



sensitivity of fluxes to all the mixing angles including CP phase 6. The signiflcant sensitivity 
to ^23 is expressed by the broad coverage in $*(z^e) and ^^{i'^) of each colored stripe for flxed 
612 and ^13. Each stripe is a bundle of curves for 6 values from to vr. The top line of the 
stripe for $*(z/e) indicates 6 = 0, while the bottom line of the stripe for ^^{i^^) indicates 
6 = 0. Next, the sensitivity to Ou was estimated only for 5 = and 6 = tt and is illustrated 
by three kinds of lines outside the colored. It is clear that the 612 does not affect the fluxes 
as much as the 6'23 does. The 613 dependency in the fluxes is illustrated by the widths of 
the flux stripes. A remarkable improvement in precision of 613 will be obtained from near 
future oscillation experiments, e.g., Daya Bay [26], Double CHOOZ [2']^, RENO [2^, Noz/A 



29|] and T2K [30|]. In FIG. [T] and FIG. [21 two values of ^13 are used for comparison: 0.20 
which is the current upper bound at 3cr CL and 0.082 which is the upper bound at 3a that 
RENO will accomplish. It is clearly illustrated that narrowing the stripes by a small bound 
of ^13 makes identifying the relative composition more realistic. 

If a series of LBL oscillations, e.g., Noz/A and T2K, are successful in achie ving the aimed 



precision level, reducing the relative range of sin^ ^23 at 3a from 79% to 42% J3l| , the range 
in $*(z/e) from the pure pion source is completely separated from the range in $*(t'e) from 
the pure muon-damped source. The estimation of composition rate in a mixed beam of two 
pure sources requires the determination of CP phase 6 as well as improved precisions of 
other parameters. 

The curves of neutrino fluxes with respect to energy show both saturation and transition. 
For instance, the initial flux ratio $°(z/^)/$°(z/e) of atmospheric neutrinos saturates to 2 



at low energy limit 9| lOj, while the expected ratio <l'*(z/^)/$*(z/e) of high energy cosmic 
neutrinos to be detected saturates to 1 at low energy limit and to 1.8 at high energy limit 
|6|]. If we can distinguish the saturation from the transition in flux - energy plots, the 
identification of neutrino beam from pure pion source or pure muon-damped source may be 
possible only after the accumulation of sufficient data of cosmic neutrino detection. 

IV. LBL: A COMPLEMENT TO TELESCOPE 

The probability P^e in Eq- © depends on cos{S + Aml-^^L/AE^), while $*(i/) in Eq. ([7]) 
depends on cos 5. If $*(z/) and P^e are the orthogonal axes as in FIG. [3|, the locus of $*(i/) 
and P^e for 6 from to 2tt completes a closed path. In FIG. [3], every locus in each figure 
passes a point in P^g — "^'(^e) space. The phrases P(t'^ -^ Ve) and ^{v^ + Va) in figures 
are equivalent to P^e and $(z/q) in text. Especially in these figures, the probability P^e 
is an expected measurement of a super beam from J-PARC to Super-Kamiokande (T2K 



collaboration). Thus, input values are adopted from Ref. 30|]. The baseline is 295km and 
the energy of neutrino is IGeV. The matter effect A in Eq. (5) is obtained with the Earth's 
density 2.8g/cm^ and the electron fraction 0.5. The Eq. (j4]) and Eq. ([5]) include only the 
leading terms. 
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FIG. 3: Eight loci for 6 G [0, 2it] passing a common point in the $*(iv) — P^e plane. 

A. Eightfold degeneracy 

A number of strategies are proposed to solve the degeneracy problems by using the results 
of LBL oscillations and reactor neutrino oscillations in future 19|-|24[. An example is to 



analyze long-baseline experiments over different oscillation distances 13], and another is to 



combine the results of reactor oscillations and the LBL |2_0']. This section checks possible 
resolutions obtainable from the combined analysis of results at telescopes and results at 
LBL. 

The leading term of P^e in Eq. ([5]) with a remarkable sensitivity to ^13 does not distinguish 
Am|^ > from Arrig^ < unless the matter effect is significant, and any detected flux $*(z/) 
is fully blind to Am|]^. However, the locus of P^e and $*(i^), as in FIG. [3]and FIG. HI reveals 
the sign of mass-squared difference Am^i due to the phase difference Am'^iL/AE^. 

Each single curve in FIG. |3] is obtained by a set of speciflc values of 623, 612 and Arn^i, 
while 6 runs from to 27r. The value of 612 is common for all curves in a flgure and each of 
623 and Arrig^ has two choices: that 623 may be a certain angle or its complementary angle 
and that Am|]^ may be |Am|;^| or — |Am|]^|. Then, two sets of {0i3,6) exist for every set of 
(6'23, Aml^), so that total eight pairs of {613, 5)'s describe the point in P^e — ^^{^e) [24]. For 
instance, the point (^'(z/e), -P/^e) = (0.33, 0.010) in FIG. [3] intersected by the eight curves can 
be specifled by the following eight pairs of (^13, 6) which belong to 'a - h'curves: 

a: (0.104,2.3) b : (0.087,3.8) 

c: (0.094,2.4) d : (0.079,3.6) 

e: (0.086,0.86) /: (0.071,5.6) 

g: (0.103,1.01) h : (0.084,5.4). 
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The above values of (6'i3, S) are divided into two groups, according to 623: [a, b, c, d] for 6*23 = 
Tc/A — e and [e, f, g, h] for 623 = 7r/4 + e where e is a positive small angle. They can be divided 
also according to the sign of Arrigj^: [a, b, e, f] for Ar7i|^ > and [c, d, g, h] for Am"^^ < 0. 
In other words, the LBL oscillation probability P^e(^23 — ^/4, sgnlAm^^), {613, 6)) = 0.010 
is degenerated by the following eight possibilities 

P^le{-(^, +, a) = P^e(-e, +, b) 

= P^e(-e, -, c) = P^e(-e, -, d) 

= P^e(+e, +, e) = P^e(+e, +, f) 

= P^e(+e, -, g) = PMe(+e, -, h), 

where P^e{~^) +, a), as an example, means that the probability is obtained when 623 — 
7r/4 = — e, sgn{/\m\^) is +, and (^13, 5) is the point on the locus 'a'. Likewise, the flux 
at telescope $*(z/e)(^23 — 7r/4, sgn{A.m\^), {613,8)) = 0.33 is degenerated by the following 
eight possibilities 

$*(z/,)(-e, +, a) = ¥{u,){-e, +, b) 

= $*(z/e)(-e, -, c) = ^\u,)i-e, -, d) 

= $*(z/e)(+e, +, e) = ^\iy,)i+e, +, f) 

= $*(z/e)(+e, -, g) = <^\ue){+e, -, h). 

Thus, a point representing data from two experiments, e.g., the point ($*(z/e), P^e) = 
(0.33,0.010) or the point (<l>*(z/^), P^e) = (0.34,0.010) in FIG. [3l always belongs to many 
curves specified by different combinations of the parameters. The number of the loci passing 
a point represents the order of degeneracy. In fact, the order of degeneracy becomes infinite 
if Arrig]^ or 623 is allowed within a continuous range. 

B. Breakable degeneracy by high energy cosmic neutrinos 

FIG. m displays the loci of P^e and $*(z/) for 612 from 0.51 to 0.69 and 6 from zero to 
27r at ^13 = 0.010 and 623 = vr/4. The projections of the two bands onto $*(z/) axis overlap 
completely with each other, and even the projections of the two bands onto P^e axis overlap 
mostly with each other. So the individual measurement, either $*(z/) or P^e, is not sensitive 
to the sign of A?7i|^. However, the combination of two measurements as in FIG. [3] or FIG. 
m is significantly sensitive to whether Amgj^ > or Amgj^ < 0. On the other hand, if 623 

11 




0,320 0.325 0.330 0.335 0.340 0.345 

<t^iVg+Vg) from pion source 



0.320 0.325 0.330 0.335 0.340 0.345 

**(>>+>>) fi'oiii pioii source 



FIG. 4: Dependence of curves on Ar7T,|^: 9i2 is allowed within 3cr CL. 

and ^13 also are allowed within certain ranges rather than with fixed values, it is hard to get 
such distinct plots according to the type of mass hierarchy. Thus, an analysis with a fixed 
6i3 as in FIG. H] can be effective only after a series of reactor or LBL neutrino oscillations, 
which will be launched ahead of neutrino detection at telescopes [26|- 30|. 

The leading term of 1 — P^^ in Eq. (jl]) is sensitive to 6*23 but does not distinguish 
^23 = 7r/4 — e from 623 = vr/4 + e. On the other hand, the flux $*(z/e) in FIG. 0, which is 
more sensitive to 623 than the probability P^e is, shows a distinction whether 623 = ir/A — e or 
^23 = 7r/4 + e. Thus, $*(z/e) which is a completely single- valued curve of 623 in the currently 
allowed range as illustrated in FIG. [1] or in FIG. [5] avoids the degeneracy caused by the 
complementarity angle of 6*23, unlike 1 — P^^ in Eq. (jl]). 

C. Source identification 

FIG. [5] is drawn under the assumption which is that the neutrino beam was purely from 
the pion source and purely from the muon-damped source. In reality, the initial flux in 
an astronomical burst is difficult to identify, whether pion is produced in py collision or pp 
collision and whether the initial condition is pion source, muon-damped source, or a mixture 
of them when p7 collision is dominant. Unlike the detection of photons at optical telescopes, 
the detection of neutrinos at telescopes should consider the change in neutrino fiavors from 
the original outbursts. Moreover, broad uncertainties at present in neutrino mixing angles 
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FIG. 5: The loci for b e [0, 27r] at 6*13 = 0.15, 0.10, or 0.05 for each B 



23- 



0,05 
-S 04 

T 

' 0,03 
.a. 

R, 0,02 
0,01 



^3=0. 97 




e23=0.62 



013=0.15,5=0 



0.20 0.25 0.30 0.35 



Am?,<0 



^3=0.97 




613=0.15,5=0 



023=062 



0.20 0.25 0.30 



).35 



*(>'e+n 



FIG. 6: The loci for ^23 G [0.62, 0.97] for various compositions. 

obstruct telescopes in describing the initial condition of neutrino beams. 

The LBL oscillation probability P^e may again share a role in specifying the source of 
cosmic neutrinos. The shaded regions of both panels in FIG. [5] show the same area which 
both sources cover commonly. So the measurement of $*(t'e) within the shadow may indicate 
a beam from a pion source (right panel) for ^23 = 0.970 or a beam from a muon-damped 
source (left panel) for ^23 = 0.623. The distinction between them may be obtained by 




fe=[)623. — fts=Ji/4. -— ft3=[1.97D 



FIG. 7: The maximum portion of pion source in a mixed beam which can be distinguished from 
the pure pion source. 
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considering P^e together with $*(t'e)- For instance, Fig. 5 tells us that P^e > 0.025 in 
the shadow can be compatible with $*(z/e) from a pion source but not with $*(z/e) from a 
muon-damped source. 

FIG. in] is an example that predicts the relative composition between a pion source and 
a muon-damped source, where 612 = 0.59, ^13 = 0.15, and 6 = 0. The change in color of 
the curves from the blue (rightmost) to the red (leftmost) indicate the change in composi- 
tion from 100% pion source to 100% muon-damped source. Unfortunately, the distinction 
of the composition by every 10% as in Figl6] lies far beyond the practicability, since cur- 
rent sensitivities of detectors are too low to use the strategy in the figure for the precise 
comparison. 

The detected flux which is partially from pion source and partially from muon-damped 
source is expressed in terms of x, the portion of pion-source flux in the total detected u^ flux 

$*(z/e) = X%{U^) + (1 - X)<l>^d(z/e), (11) 

where $p is a flux from pion source and ^md is a flux from muon-damped source. Even 
though the exact value of ^13 was known, the undetermined 6 limits the composition of 
mixed flux or the pure pion-source flux to have blind ranges. FIG. [7] describes the maximum 
portion X of $p(t'e) in a mixed $*(t'e) for various ^13 and 623 that can be distinguished from 
the pure $p(z/e), where < x < X. The reason of decreasing X as 6^13 increases is because 
amplitudes of 6 curves become more sizable. The upper(blue) lines are affected only by the 
range in 6, while the lower(red) lines are affected by systematic uncertainties as well as by 
6. The sensitivity of IceCube detector to astrophysical source was discussed for neutrinos at 



TeV to PeV energies 32] ■ When ^^ is proportional to E ^, the systematic uncertainty is 



+10/ — 15%, and when ^^ is proportional to E ^, the systematic uncertainty is +5/ — 20%. 



-2 



In case of E~ spectrum, the composition of a distinguishable mixture from the pure pion- 
source beam is drawn by the lower (red) curves in FIG. [61 Within the current sensitivity 
of the IceCube detector, the mixed flux with more than 40% (or 20%) pion source cannot 
be distinguished from the pure pion source when 623 is 0.970 (or vr/4) and 613 reaches the 
3cr upper bound at RENO. Within the current upper bound of 6^13 as seen in the flgure, the 
pure pion source cannot be distinguished even from the pure muon-damped source, unless 
P^e does help the identification in such a way as in FIG. [5] and FIG. El 
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V. CONCLUDING REMARKS 

Neutrino telescopes will detect neutrino beams out of astronomical bursts. However, the 
results to obtain are barely helpful in describing the initial condition of cosmic neutrino 
beams since the uncertainties in neutrino masses and mixing angles are broad and the fluxes 
to be measured are sensitive to the masses and the mixing angles. On the other hand, the 
improvement of the precision in parameters is also hard to attain by using results at neutrino 
telescope itself since even an original beam as initial condition cannot be defined without 
using the telescope. 

We took a strategy to consider the fluxes to be detected at a telescope like IceCube in 
company with oscillation probabilities at a LBL T2K. The expected fluxes are examined 
for the sensitivities to mixing angles, in comparison with the sensitivities of other types of 
oscillation probabilities to mixing angles. A few restricted cases were presented as examples 
to show that neutrino fluxes at telescopes may be useful to resolve the degeneracies embedded 
in terrestrial neutrino oscillations. It was followed by the discussion on the limit of source 
identification which is allowed within the sensitivity of IceCube to astrophysical neutrinos. 
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